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Abstract. The isospectral flows of an n th order linear scalar differential opera- 
tor L under the hypothesis that it possess a Baker-Akhiezer function were originally 
investigated by Segal and Wilson from the point of view of infinite dimensional Grass- 
manians, and the reduction of the KP hierarchy to the GePfand-Dikii hierarchy. The 
associated first order systems and their formal asymptotic solutions have a rich Lie al- 
gebraic structure which was investigated by Drinfeld and Sokolov. We investigate the 
matrix Riemann-Hilbert factorizations for these systems, and show that different fac- 
torizations lead respectively to the potential, modified, and ordinary GePfand-Dikii 
flows. Lie algebra decompositions (the Adler-Kostant-Symes method) are obtained 
for the modified and potential flows. For n > 3 the appropriate factorization for 
the Gel'fand-Dikii flows is not a group factorization, as would be expected; yet a 
modification of the dressing method still works. A direct proof, based on a Fredholm 
determinant associated with the factorization problem, is given that the potentials 
are meromorphic in x and in the time variables. Potentials with Baker-Akhiezer 
functions include the multisoliton and rational solutions, as well as potentials in the 
scattering class with compactly supported scattering data. The latter are dense in 
the scattering class. 



o 



1. The Gel'fand-Dikii hierarchies 



Gel'fand and Dikii [GD] constructed a hierarchy of isospectral flows of the n th 



order scalar differential operator 



3=0 



u 1 (x)D n - J ; D = -4~; u = l, ui=0. 

% ax 



where uj = Uj(x), j > 2, are elements of the Schwartz class «S(R). The flows are 
given by 

L=[L k J n ,L], k=l,2,... (1.1) 

where fc ^ mod n, and L k ^ n denotes the differential part of L k l n considered 
as a pseudodifferential operator. We shall refer to the coefficients Uj of L as the 

potential. The coefficients of L k ^ n are obtained by solving formal recursion relations 
and are universal differential polynomials in the {«?}• 

The forward and inverse scattering problem for L has been investigated for po- 
tentials in iS, and the Gel'fand-Dikii flows for potentials in this class have been 
constructed by Beals [B], and Beals, Deift, and Tomei [BDT]. Action-angle vari- 
ables for flows in S were constructed in terms of the scattering data by Beals and 
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Sattinger [BS1]. Segal and Wilson [SW] investigated these flows under the ansatz 
that L possess a wave function of Baker-Akhiezer type: w(x, z) exp{ixz}, where 
w is analytic at z = oo and takes the value 1 there. From now on, we denote 
this class of potentials by B. They analysed the flows from the point of view of 
infinite dimensional Grassmanians, and the reduction of the KP hierarchy to the 
Gel'fand-Dikii hierarchy. 

The n th order scalar equations can be converted tonxn first order systems. 
These have a rich Lie algebraic structure (cf. Drinfeld and Sokolov [DS]), just as 
the Dirac equation rather than the Klein Gordon equation exhibits the internal 
symmetries of the electron. All reductions are gauge equivalent, but besides the 
Gel'fand-Dikki flows, one finds others, such as the modified and potential Gel'fand- 
Dikii flows. We analyze the matrix Riemann-Hilbert problems corresponding to the 
first order systems. For potentials in £>, the formal asymptotic series in [DS] are 
actually convergent. 

Given a Baker-Akhiezer function, it is easily seen that L has a basis of such wave 
functions. On the other hand, there is always a basis of wave functions which are 
entire in z. Denoting the Wronskian of this basis by <j> + and that of the basis of 
Baker-Akhiezer functions by 0_ we form the quotient g = <f)Z 1 4>+- The matrices 
4>± can be chosen in such a way that g is a function only of A = z n . Conversely, 
appropriate Riemann-Hilbert factorizations of a given a matrix g(X) together with a 
standard argument (the so-called "dressing method"), yields a potential in B. Thus, 
while we cannot precisely characterize the class of potentials £>, we nevertheless can 
construct them by solving a factorization problem. 

All conversions of the n th order equation to a first order system are gauge equiv- 
alent, and the entries of L are differential polynomials in the entries of the potential 
q [DS]. The canonical conversion (by Wronskians) leads to the Gel'fand-Dikii flows 
themselves. The modified Gel'fand-Dikii flows are obtained by first factoring L as 
a product of n first order operators, and constructing a system of AKNS type. We 
investigate in particular the standard and modified Gel'fand-Dikii flows together 
with a third type we call the "potential" flows. All three are obtained by factoring 
exactly the same matrix g, but taking the factors in different infinite dimensional 
submanifolds. For the potential and modified flows, that submanifold is a subgroup 
and the components of the connections are obtained by Lie (or loop) algebra de- 
compositions ( a la Adler-Kostant-Symes) . For the Gel'fand-Dikii flows themselves 
the submanifold is not a group for n > 3. A Lie algebra decomposition is not 
available, but the dressing argument need only be modified slightly. 

Schilling [Sch] obtained a Lie algebra decomposition for the first order system of 

equations for <f>j = v, j = 1, . . . , n, but did not consider the corresponding 

factorization problem and flows. 

The hypothesis that L possess a Baker-Akhiezer function is a very strong one, as 
observed in [S W] . If q is rational in x the wave functions are in general multiple val- 
ued in the punctured neighborhoods of the poles of q; but if L has a Baker-Akhiezer 
function then the wave functions have trivial monodromy (see §2.3). Moreover, if 
q is a priori defined only on some interval I on the real line, and L has a Baker- 
Akhiezer function for x G /, then q has a meromorphic extension to the entire 
complex plane (Corollary 3.1.7). 

In the case n = 2, assuming, say, that (1 + |cc|)u2 G the assumption implies 

j-U„j- j-i a — — £E„; — + :„ — ™„„„+u. „ +„j t„ ;f :„ — „i 
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in z, then the reflection coefficient vanishes, and we have a so-called "reflectionless 
potential." Properly interpreted, the same holds for general n, so that in a sense the 
ansatz is a generalization of the notion of reflectionless potential. Segal and Wilson 
showed that the class B includes the multi-soliton solutions, the rational solutions, 
and in general the "algebro-geometric" solutions of the KdV equation. We show in 
this paper that the class B is actually somewhat larger, and in fact that Sq D B is 
dense in Sq (Theorem 2.1.7). (So is the class of potentials for which the scattering 
transform is well-defined in [BDT]). 

We show, by arguments based on the Fredholm determinant of an operator of 
trace class associated with the factorization problem, that the potentials in B are 
meromorphic in x and in the time variables under the Gel'fand-Dikii flows. Such 
results were proved by Segal and Wilson [SW] using the theory of the r function 
for the Kadomtsev-Petviashvili (KP) hierarchy, together with the reduction of the 
KP hierarchy to the Gel'fand-Dikii flows. Our proof is similar in spirit but requires 
considerably less machinery; furthermore, it can be extended to other systems in 
the Drinfeld-Sokolov hierarchy, as well as the n x n AKNS flows. 

The non- vanishing of the Fredholm determinant for the system guarantees both a 
left and right factorization; whereas the non-vanishing of the r function guarantees 
only a one-sided factorization. 

The factorization technique used here provides an efficient computational algo- 
rithm for constructing the rational and multi-soliton solutions. We hope to take 
this up in a future paper. 

The restricted class B is distinct from the potentials for which the scattering 
transform of L is well-defined. The "scattering class" of potentials is invariant 
under perturbations by functions in the Schwartz class S, whereas B is not. On the 
other hand, B includes the rational solutions, which, due to their slow decay, or to 
possible poles on the real axis, are not in the scattering class. 

For generic potentials in the scattering class, the wave functions are meromorphic 
in sectors of the complex plane of angle rr/n. The scattering class leads to a 
Riemann-Hilbert factorization problem on a set of rays S running from the origin 
to infinity, [BDT] while potentials in B are described by a factorization problem on 
a circle centered at the origin. This factorization problem is considerably simpler 
than that associated with the inverse scattering problem; moreover, it is closely 
tied to the infinite dimensional Lie Group/ Lie algebra structure associated with 
the first order systems [DS]. 

In §2 we develop the forward problem for potentials in B by showing that the 
ansatz leads to a matrix factorization problem on a circle centered at the origin in 
the complex A plane. 

In §3 we consider the inverse problem, formulated as a Riemann-Hilbert problem. 
The parameters x and t appear analytically in the factorization problem, and the 
solutions of the factorization problem are shown to be meromorphic in the parame- 
ters. The dressing method is used to construct a flat connection from the solutions 
of the factorization problem and show that the resulting connection potentials are 
meromorphic in x and t. Similar considerations arise in the proof of the Painleve 
property for solutions of isomonodromy deformation problems [Mai] (cf. also Miwa, 
who proved the Painleve property using the r function.) 

In §4 we discuss the infinite dimensional Lie structure of the factorization prob- 
lem, using gradings of the algebra introduced by Drinfeld and Sokolov. In §5 we 
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Gel'fand-Dikii flows. 

This work was carried out during the second author's visit to the IMA and the 
School of Mathematics in the summer of 1991. He would like to thank the IMA and 
the University of Minnesota for their hospitality. Both authors acknowledge valu- 
able discussions with R. Beals, J. Dorfmeister, and I. Herbst concerning various 
aspects of scattering theory, factorization theory, and the Gel'fand-Dikii hierar- 
chies. A number of techniques used in this paper came out of joint work between 
the first author and Richard Beals on the Gel'fand-Dikii flows and isomonodromy 
deformations. 

We also benefitted from comments by the referee and H. Flaschka. 



2. The forward problem 

2.1 Reduction to a first order system. The scattering theory of L proceeds 
by converting the n th order scalar equation Lf = Xf to a first order system. Let 
4>j = Di~ l f for j = 1, . . . , n. Then the column vector <fi = (0i, 2) • • • ; <PnV satisfies 
the first order system 

(D — J\ — q)cf) = 0, (2.1.1) 

where 



Jx 



/0 1 

1 



Va o o 



\ 



1 

J 







Q = - 












U n -1 



\ 





U 2 / 



The linear space of all such q constitutes a special class of potentials, which 
we call the Gel'fand-Dikii potentials. For the rest of this section, q will denote a 
potential in this class. In general, we do not assume the entries of q are in the 
Schwartz class; but we say that q G «S(M) if each entry of q belongs to S. 

Matrix solutions of (2.1.1) are necessarily Wronskians, the entries of the first row 
being scalar solutions of Lf = Xf. 

Let f(x,z) be a Baker- Akhiezer function for L, dj, 1 < j < n, the n th roots 
of unity, and let fj(x,z) = f(x,ajz). The Wronskian matrix W = \\D fj\\, 
1 < k, j < n, satisfies (2.1.1). From the asymptotic behavior of / we conclude 



W H|£> fc ~V XQ ^|| = A z e ixzJa = d(z)e lxzJ ° 
z ^oo 



A a 



where 



and 



? a n)i 


d(z) 


= diag(l,2, . . 






( 1 


1 1 ... 


1 








a n 


A a = 


a\ 




a 2 



\ 



A z = d(z)A c 
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The matrix function ip = d l W K a l satisfies 

L>V = (zJ + q(z))i(> lim ^e~ ixzJ = 1 (2.1.2) 

z^oo 

where q(z) = d~ 1 qd, d~ x J\d = zJ, and 



J = 



/0 1 ... ^ 

1 ... 

... 1 

Vl ... J 



This proves: 

2.1.3 Theorem. L has a Baker- Akhiezer wave function iff (2.1.2) has a funda- 
mental solution of the form ip = m(x, z)e %xzJ where m(x,z) is analytic at z = oo 
and tends to I there. 

Since tr(zJ + q) = it follows that detip, and hence detm is constant in x. Since 
m — > I as z — > oo, detm = 1 at infinity. We may therefore renormalize it so that 
det m is identically 1 in a neighborhood of infinity. From now on we shall refer to 
m as a Baker-Akhiezer function for (2.1.2). 

We now discuss, briefly, the scattering data of L when q e S [BDT] . Let E = {£ : 
Re ix£(otk — dj) = 0, j, k = 1, ... , n}, and denote the 2n sectors which complement 
£ by O^, = 1,2,..., 2n. Solutions of (2.1.2) are constructed in each Vt v for which 
the corresponding m v satisfies 

Dm = [zJ, m] + q(z)m, (2-1.4) 

and the "radiation conditions" 

lim m(x,z)=I, sup \m(x, z) \ < +oo, (2.1.5) 

x^ — oo x 

where | • | denotes the sum of the absolute values of the entries of the matrix m. The 
wave function m v is uniquely determined, except possibly at a bounded discrete 
set Z, where it has poles. The points of Z may cluster only along the rays of E. 
Away from those cluster points, m v is continuous up to the boundary of Vt v . For 
a generic set of potentials, denoted by Sq, Z consists of a finite set of points. For 
these q the boundary values of m exist everywhere on E\{0}. 

The scattering data consists of the multiplicative jumps (continuous component) 
of m across the rays of E, together with its principal part (discrete component) at 
points in Z. The continuous component of the scattering data is algebraically 
related, in the case n = 2, to the reflection coefficient. The multi-soliton potentials 
of L = D 2 +u are examples of potentials in BdSo with a trivial multiplicative jump. 
These are typically referred to as "reflectionless potentials". We show below that 
for potentials in B fl So the multiplicative jumps are compactly supported. Thus B 
is an extension of the class of potentials with compactly supported scattering data, 
and in particular contains the reflectionless potentials. 

Let J denote the kernel of ad J, and let O denote a neighborhood of infinity in 
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2.1.6 Lemma. Let O be a neighborhood of infinity in the complex plane, let mi 
andmi be two solutions of (2.1.4) which are analytic in OC\Vl v , bounded for x G K 
and, for all real x, tend to I as z — > oo in Q, v . Then mi = rri20~v{z), where 5 V 
takes its values in J , is analytic in O D Cl u , and tends to I as z tends to infinity 
along rays in the interior ofVl v . 

If mi and m2 satisfy (2.1.4) for z G O, and take the value I at infinity, then 
mi = m2S(z) in O where S(z) takes its values in J , is analytic in O , and takes the 
value I at infinity. 

Proof. The matrix w = mj~ 1 m2 satisfies Dw = [zJ,w], hence w = e lxzJ 5{z)e~ lxzJ 
for some matrix valued function 6 which is independent of x and analytic in OC\VL v . 
Since w is bounded for — oo < x < oo and tends to I as z tends to infinity, 5 must 
have the required properties. This proves the first statement. 

In the second case, w is analytic at infinity, so has a Laurent expansion T,jWj(x)z~- 
where the coefficients Wj satisfy the recursion relations [J, Wj+i] = dwj/dx, with 
wq = I. It follows that [J, wi] = 0, so that w\ G J . Since ad J is semi-simple, its 
kernel and range intersect only at zero. Hence [J, W2] = dwi/dx and w\ G J imply 
that wi is independent of x and W2 G J. By induction, all the coefficients Wj are 
independent of x and lie in J. □ 

We are now ready to prove 

2.1.7 Theorem. In the following, let q G S. 

i. If the continuous component of the scattering data is trivial, and Z is finite, 
then L has a Baker- Akhiezer function for which m is rational in z and is 
bounded for all real x for regular values of z. 

ii. Conversely if L has a Baker Akhiezer function for which m is rational in 
z and bounded for all real x for regular values of z, then the continuous 
component of the scattering data is trivial. 

hi. If the scattering data of q is compactly supported, then q G B; and con- 
versely, if L has a Baker- Akhiezer function for which m is bounded in x for 
regular values of z then the scattering data is compactly supported. 

iv. So n B is dense in So- 

Proof. To prove i., define a sectionally meromorphic function m by m(x,z) = 
m l/ (x,z) for z G O^, where m v are the wave functions constructed in the for- 
ward scattering problem. Since m has no jumps across £\0 and Z is a finite set 
m is meromorphic in C\0. (It takes the value / at infinity.) By Theorem 8.15 in 
[BDT], z n ~ 1 d(z)m(x, z)d~ 1 (z) has a finite limit at z = 0. Consequently, m is a 
rational function of z which is bounded in x for regular values of z. 

Conversely, let m and m v denote the fundamental solutions of (2.1.4) obtained 
respectively from the Baker- Akhiezer function and the forward scattering problem. 
By the first statement in Lemma 2.1.6, m = m v 8 v {z) for zGfi„, where b~ v takes its 
values in J and is analytic in the intersection of Vt v with a neighborhood of infinity. 
Since m v is continuous up to the boundary of Q u , 

m~ l m v+ i{x,i) = SvS'l^Z), £ G E v , 

where = n Q, v +\. Since ad J is a derivation, s„(£) = SvS„+i G J. 

We wish to show that each s u is the identity matrix. We can just as well work 
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denote the roots of unity. In that case, s v is also diagonal, and the above equation 
may be written as 

m v+1 (x,£) = m„(x,£)s u (£), 
or, since s v commutes with J, 

e ix * J m v+1 e- ix * J = e ix ^ ' m v e~ ix ^ ' s v {i) (2.1.8) 

For £ G £, let II^ denote the projection 



ejk if Re i£{ctj — a k ) = 
otherwise, 



where ej k denotes the matrix with a 1 in the jk place and zeroes elsewhere. We 
prove below that 

lim U,e lx ^ J m K e- ix ^ J = I + v K , k = v,v + 1, (2.1.9) 

x^ — oo 

where v v , v v +\ lie in the range of 11^ and are respectively strictly upper and lower 
triangular in an appropriate ordering of the roots otj. From the definition of II^, we 
see that U^as = (U^a)s for a diagonal matrix s. Applying II^ to (2.1.8) and taking 
the limit as x — > — oo, we obtain 

I + V v+ i = (I + v v )s v . 

By the structure of v u and v u +i we must have s„ = I and v v = v„+i = 0. 

We establish (2.1.9) by considering an integral equation for the boundary values 
of the wave functions. For z G 0„\Z, m v satisfies the Fredholm integral equation 
[B] 

m v (x, z) = I+ [ (nj + IUy)e iz ( x -rt J (qm v )e- iz ( x -rt J dy 

J — oo 
poo 

- / Tlte lz{x - y)J (qm u )e- lz{x - y)J dy, 

J X 

where Ilo projects onto the diagonal matrices and IlJ denote the projections 



°jk 



ejk if ± Re iz(aj — a k ) > 0, 
if ± Re iz(dj — a k ) < 0. 



An ordering of the otj may be chosen so that, for z G 0„, II f are the projections 
U± onto the strictly upper and lower triangular matrices. Letting z tend to £ G H v 
from 0„, we obtain the following integral equation for the boundary values of m v : 

m v (x,S) = I+ [ X (U-+U )e i ^ x -^ J (qm^e- i ^ x -^ J dy 

J -°° (2.1.10) 

"OO 



/•OO 

/ U^ x - y)J '(qm v )e-^ x - v)J ' dy. 
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Conjugating (2.1.10) by e lx ? J ? applying 11^, and taking the limit as x — > — oo, we 
obtain 

/oo 
U,U + e~^ yJ (qm^ yJ dy = I + v v . 
-oo 

For 2 G 0„+i, 11+ no longer projects onto the upper triangular matrices; however, 
11^11+ = U^U-. Hence, repeating the previous arguments, we obtain 

/oo 
UM-e-^ yJ (q mi , +1 )e^ yJ dy = 1 + v v+1 . 
-oo 

This establishes (2.1.9), hence ii. 

If the scattering data is trivial outside some compact set, then there is some 
open neighborhood O of infinity on which m u = m u+ i in Y* v n O, and we obtain a 
Baker-Akhiezer function by setting m = m u for z G 0„ n O. Hence q G B, and iii. 
is established. 

The result in iv) is a consequence of the fact that potentials with compactly 
supported scattering data are dense in So in the Schwartz topology. This fact is 
stated without proof in [BDT] (cf. p. 154). A rough sketch of the proof runs like 
this. Let q G So, an d denote its scattering data by v. (We also include the discrete 
data in v.) On E each entry of v — I belongs to <S(E) ([BDT], Theorem 13.1). In 
particular, each entry of v — I is rapidly decreasing at infinty and has an asymptotic 
expansion in positive powers of z as z — > along each ray E^ of E. We approximate 
v by smooth data v n such that v n = v in \z\ < n, and v n — I vanishes identically 
in the region \z\ > n + 1. Since the inverse transform is defined in a sufficiently 
small neighborhood of v, it is defined for each v n for sufficiently large n; and the 
corresponding potentials q n belong to B. The inverse transform, as constructed in 
[BDT] , is the composition of two continuous processes: a rational approximation of 
the scattering data; and the inversion of a Fredholm operator close to the identity. 
Hence the inverse transform is (locally) continous from Z/2(E) to Z^Q^)- Using the 
estimates in [BDT], one could also prove, with some effort, that the corresponding 
potentials q n converge to q in the topology of S. □ 

For potentials in S the m v are uniquely determined by their normalization as 
x — > — oo. No such unique normalization is possible for the Baker-Akhiezer wave 
functions; but we do obtain: 

2.1.11 Theorem. If q G B then a fundamental solution of (2.1.2) can be con- 
structed which satisfies the symmetry condition ip{x, az) = d~ 1 (a)ifj(x, z)d(a), whre 
a is a primitive n th root of unity. 

Remark. This symmetry condition does not uniquely determine the wave function. 

Proof. First note that (2.1.4) is invariant under right multiplication by holomorphic 
S(z) G J which tend to I as z tends to oo. These matrices have the representation 

n-l 

8{z) = 5j(z) J\ &j{ z ) ~ ¥ $jo as z — > oo. 
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Since d~ l {a)zJd(a) = azJ and d~ 1 (a)q(z)d(a) = q(az), (2.1.4) is also satisfied 
by d{a)m{x, az)d~ 1 (a). By the second statement of Lemma 2.1.6 there is a 8{z) 
with values in J such that 

d{a)m{x, az)d(a)~ 1 = m(x, z)S(z). 

Let q be regular at some point xo, and let a(z) = m(xo, z). Then 

a(az) = d- 1 (a)a(z)5(z)d(a). (2.1.12) 

The function m is uniquely determined by its values at xq. Moreover, m/i is also a 
solution of (2.1.4) for any holomorphic matrix /j,(z) with values in J. For the proof 
of Theorem 2.1.11 it therefore suffices to find a matrix n(z) with values in J such 
that 

d(a)a(az) n(az)d~ 1 (a) = a(z)/j(z); 

for it then follows that m(x, z)[i satisfies the required symmetry condition. Usi;ng 
(2.1.12) we may write this condition as 

8{z) = /j(z)d(a)^- 1 (az)d- 1 (a) (2.1.13) 

where 

8, /ie J, 8(oo)=I, det<5 = l. 
We represent 8 and // in the form 

n— 1 n— 1 

8{z) =exp{5> j (*)J'}, n(z) = exp{^(^}. 
i=i i=i 

Iterating (2.1.13) n times, we obtain 

/ = 8(z)d(a)8(az) . . .8(a n - 1 z)d(a). 

Bringing each of the n matrices d(a) all the way to the left in the above expression 
(and noting that Jd(a) = d(a)aJ), we find 

n— 1 n— 1 

^exp^^^a-V'}, 

j = l k=0 

and therefore 

n-l 



J2vj(a k z)a k = 0, j = l,...,n-l 



A:=0 



These equations represent a set of constraints on the r]j : which are given functions. 
From (2.1.13) we find by similar arguments, 

Vj (z) = 9 j (z)-a-i9 j (az). (2.1.14) 

We seek a solution of (2.1.14) in the form 

n-2 

d i( z ) = ^2 v kVj(® h z), 3 = l,-..,n-l, 
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where the coefficients Vk (more precisely Vj^) are to be determined. The three 
preceeding equations reduce to the linear system 



/ 


1 







— a _J 


1 







—a 


V 








a 



-j 




1 







a-O'+i) 



vi 



\V n -2 / 







Vo/ 



. . . -a~ j 1 + a-W" 2 ) / 
The coefficient matrix on the left is invertible, as can be seen by reducing it to 
row echelon form. Adding the first row to the second multiplied by oP , etc. one 
arrives at a matrix in row echelon form whose diagonal is 

1, a j , a 2j , - a~ j ). 

The last entry is non-zero, since a is a primitive n th root of unity. □ 

2.2 Properties of the wave function. 

We next analyze the structure of the wave function m constructed in theorem 
2.1.11. From the symmetry condition 

ip(x, az) = d~ 1 (a)ip(x, z)d(a) 



and the expansion 

we deduce 
We claim that 



rn 



j=0 



m 



d (a)mjd(a) = a J mj, j > 0. 



(2.2.1) 
(2.2.2) 

(2.2.3) 
(2.2.4) 



rrij(x) = Ej(x)J J , Ej diagonal, E = I. 

First, from the identity = eaJi~ % we see that any A e M n (C), the space of 
n x n matrices with complex entries, can be represented as 

n-l 

A = J2 E kJ h , E k diagonal. (2.2.5) 

Furthermore, 

<r 1 (a)J k d{a) = {aJ) k , k e Z. (2.2.6) 

Expanding rrij in the form (2.2.5) and using (2.2.6) we see that Ek = for j ^ 
k mod n. This establishes (2.2.4). 

From now on we often suppress the x dependence of the Ej. Further restrictions 
are placed on them by the requirement that m satisfy (2.1.4). First we note that 

n-l 

q(z) = Q j+ i(zJ)- j , Qj = diag (0, 0, ... , - Uj ). (2.2.7) 
i=i 

Substituting (2.2.2) into (2.1.4) we obtain a sequence of recursion relations for the 
rrij = EjJ~K Straightforward calculations lead to 

(a - I)E 1 =0, (2.2.8) 

j 

(a - I)E j+1 =DE J - Y,Qk+i<r- k (Ej-k), j > 1 

fc=i 

where we set Qk = for k > n and cx(Ej) = JE 3 J~ X . We shall need the following 
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2.2.9 Lemma. Let V denote the algebra of diagonal n x n matrices and let V be 
the subspace of traceless such matrices. Then 

i) a E Aut V 

ii) a n = I 

hi) Range (a — I) = Vq, Ker(a — I) = C I 

iv) On Vq and for n > 2, 

(a - I)' 1 = -(a 2 + 2a 3 + • • • + (n - 2)a n ~ 1 - I). 

Proof. Items ( i— iii) are immediate from the definition of a, while (iv) follows from 
the identity 

J + (t + (t 2 H + (T n " 1 =0 onV . 

2.2.10 Theorem. Let D\ denote the matrix diag (1, 1, . . . , 1, 1 — n). Then 

Eo =1, Ei = evil, 
j-i 

Ej =ivjl + J2(D k Uj-k)(o- - I)~ k D u j > 2. 
k=i 

The ujj are determined from the coefficients Uj of L by a sequence of recursion 
relations 

Uj = —nDujj-i + hj 2 < j < n 

where hi = and hj is a polynomial in u±, and their derivatives for 

J>2. 

Proof. By (2.2.8) E\ is a scalar multiple of the identity, so E\ = By the second 
relation in (2.2.8) 

(a - I)E 2 = DEi - Q 2 I. 

This equation is solvable provided that tr (DE\ — Q 2 ) = nDu\ + u 2 = 0; in that 

case 

DEi - Q 2 = DujiI - nDu)ie nn = Duji(I - ne nn ) = (D< j ji)D 1 , 

and 

E 2 = oo 2 I + (DuJi)(a - 7) _1 Di. 

Assume the result holds for j > 2. Since (a — I) -1 maps V to itself, trDEj = 
txDujjl = nDujj, hence the recursion relation (2.2.8) for Ej +1 is solvable provided 

j 

nDuj-tr Y,Qk+iv~ k (Ej-k) =0. (2.2.11) 

fc=i 

It follows from (2.2.11) that 



3 3 

^Q k+1 a~ k {Ej- k ) = e nn tr^2 Q k+1 a~ k (Ej- k ) = n(Duj)e nn , 
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so that (2.2.8) may be written as 

J-l 3 

(a - I)E j+1 =DujI + D k+1 u^ k {a - I)- k D x - ^ Q k+1 a- k {E 3 . k ) 

fc=l k=l 

= (Du J )D 1 +J2D k+1 0J j - k (a - I)~ k D 1 

k=i 

j 

= J2D k u J+1 . k (a-I)- k+1 D 1 . 

k=i 

Therefore 



E j+1 = + ^D k uj j+1 - k {a - I) k DL 

k=i 

For 2 < j < n - 1 (2.2.11) implies 

3 

nDojj + ^2 u k+i Tr[e nn( 7- fc (^_ fe )] 
k=i 

3-1 

= nDujj + u j+1 + ^ u k+1 TT[e nn a~ k (Ej_ k )] = 0. 
fe=i 

This yields the appropriate recursion relation for u>j and the induction step is com- 
pleted. □ 

Now we are ready to state the main theorem of this section. 
2.2.12 Theorem. Let rh = d(z)md- 1 (z). Then 
i) 

oo 



3=0 



where the Ej 's are given in Theorem 2.2.10; 

ii) rh(x, az) = rh(x, z), hence rh = rh(x, A); 

iii) rh is regular at A = oo and 



lim rh(x, A) = I + £, 



where £ is a strictly lower triangular matrix. 
iv) det rh = 1 



Proof. Observe that d(z)zJd 1 (z) = J\ and Ej diagonal imply that d(z)Ej(zJ) j d 
EjJ^ J . This proves i) and ii) follows immediately. By the definition 



n-l 
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Note however that 



lim J 7 = 



1 
Vo 











0/ 



K. 



Thus 



n-l 



m(x, oo) = E j RJ =I + Z 



where £ is strictly lower triangular. Finally detm = detm = 1. □ 

2.3 Factorization for the Gelfand-Dikii flows.. 

The existence of a Baker- Akhiezer function for L leads to a factorization problem 
on a circle centered at the origin in the A plane; specifically: 

2.3.1 Theorem. Let q e B. Then there exist wave functions <j)±(x,\) of (2.1.1) 



such that (f) + is entire and 



me tx x , where m(x, oo) = I + £, £ strictly lower 



triangular. <p± satisfy the factorization problem 

g(X) = (fiZ 1 ^, \)4> + (x, A), 



(2.3.2) 



where g is analytic in p < |A| < +oo for some p > 0. By a suitable rescaling of the 
potentials, we may assume that p < 1. 



Proof. The wave function 

4>- = d(z)tp-d~ 1 (z) = d(z)m-(x, z)d~ 



l e ixJx 



rh(x, A)e 



ixj\ 



satisfies (2.1.1) and, by Theorem 2.2.12, has all the required properties. Since 
(2.1.1) depends analytically on A, there is a unique solution <j>+(x, A), which is 
entire in A and satisfies the initial value condition 0+(O, A) = /. (We assume, for 
convenience, that q is regular at x = 0.) Since 4>± both satisfy (2.1.1), the matrix 
(pZ 1 (p+ is independent of x. 

Given a Baker- Akhiezer function for the potentials Uj(x) which is analytic in 
the exterior of some circle in the A-plane, we define rescaled potentials by u* (x) = 
s^Uj{sx), s > 0. If q s (x) denotes the matrix of rescaled potentials, then <p s = 
d(s)4>(sx, s~ n \)d~ 1 (s) satisfies (D — J\ — q s (x))<p s = 0. Replacing the original 
wave functions by the rescaled ones, we obtain a factorization problem on any 
circle in the A plane we choose. In particular, we may assume (2.3.2) holds for 
some p < 1. □ 

The isospectral flows of L are compatibility conditions for the pair of scalar 
equations 

Lf = Xf D t f = (L k J n )f, D t = -id/dt. 

By (2.1.1), all ^-derivatives of / can be written as linear combinations of the first 
n — 1 derivatives of / with coefficients which are differential polynomials in the 
Uj. The equation for the time evolution of / is therefore equivalent to a first 



order system of the form (D t — v)W 

J J ;„ Co i o; j-U„ CR„ 



0, where W = 



\D k 1 fj\ \ is the Wronskian 



j;£F + 



,1 



14 



D.H. SATTINGER & J.S. SZMIGIELSKI 



the coefficients of L, v aiso possesses such a structure. When q = we find that 
v = so the Gel'fand-Dikii flows are zero curvature conditions for the connection 
D — J\ — q, Df — J^ — Bf~, n , Bk, n denoting a matrix valued differential polynomial in 
uj of degree k — 1 in A. Let q(x, t) e B evolve according to the (n, k) Gel'fand-Dikii 
flow. The wave functions <f>± of Theorem 2.3.1 satisfy 

(D - J x - q)cp± = 0, (D t - J k x - £fc, n )0± = 0. 
They may be written as 

<j>± = m±(x,t,\)e l{xJ * +tJ *\ 

where m + is an entire function of A, and m_ satisfies the conditions of Theorem 
2.1.12. From Theorem 2.3.1 we have: 

2.3.3 Corollary. The functions rh± satisfy the factorization problem 

el{xJx+ tj*) gWe -r( x j x +tj*) = (2 _3_ 4) 

where again g is analytic in an exterior domain containing the unit circle, except 
possibly at A = oo. 

Later, in Theorem 4.3.4, we shall see that the factorization given in (2.3.4) is 
unique. Using this, we may prove: 

2.3.4 Theorem. Suppose q is rational in x in a domain U and in addition is a 
Baker- Akhiezer potential. Then the monodromy of the wave functions is trivial. 

Proof. Suppose 

e xJx ge -xJ x — m - 1 m+ 

and suppose rh± denote the analytic continuations of m± around a pole of q. Since 
the process of analytic continuation is continuous, we must have 

e xJx ge -xJ x =rhZ 1 rh + . 

By the uniqueness of the decomposition (Theorem 4.3.4 below), we must have 
rh± = m±. □ 

3. The Inverse Problem 

3.1 Factorization on the circle. 

In §2 we derived the factorization problem (2.3.4) on the unit circle, beginning 
with L and its isospectral flows for potentials in B. The inverse problem consists 
of solving a factorization problem and constructing the flat connection from the 
factors rh±. We begin with a brief summary of the solution of matrix factorization 
problems on the unit circle, (cf. [GK].) 

Let G denote the loop group 

/~i r„. .^x^i r- n/r /rm j„+ i „ :„ „^l\l ^ 
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together with its subgroups (under standard matrix multiplication) 

G + ={ g eG:g = J29kX k h 

k>0 

G-={geG:g = J2 9kX k , g = 1} 

k<0 

We always assume p < 1, so that the domain of holomorphy of g contains the unit 
circle. In practice, we use g to signify either the function g(X) or the sequence of 
Fourier coefficients {gk}- For g e G 

J2\9k\s h < +oo 

k 

for all s > p, where \g^\ denotes the sum of the absolute values of the entries of gk- 
Therefore for any s > p, G is a subset of the Hilbert space H s defined by 

H' = {g : | \g\ | 2 = ^ fj tr [<,(«,*%(«,*)*)] d6 = i- £ , 2fc tr ^ < +oo}. 
We define 



= g # s : = for k < 0}, H s _ = {g e H s : g k = for k > 0}. 

Then are the Hilbert spaces of boundary values of matrix- valued functions 
holomorphic in the interior and exterior of the unit circle. From now on we some- 
times drop the superscript s. Moreover, if cf>_ e H 8 _ for all s > p and det0_ = 1 
then <p_ e G-. 

We associate a Fredholm determinant with the factorization problem. (Segal and 
Wilson obtain the r function as a determinant associated with a Riemann-Hilbert 
problem.) The Fredholm determinant is defined for operators of the form I + A 
where A is of trace class on a Hilbert space H [Si] . The operators of trace class on 
H form a complex Banach space 3± with norm 

||A||i = trV AA*. 

The nonlinear functional det (I + A) is well defined and differentiable on Ji. Under 
the identification of 3\ with the space of bounded linear operators on H, its Frechet 
derivative is 

(I + A)' 1 det(I + A). 

The operator 

D(A) = -A(I + A)' 1 det(I + A) 
is the first Fredholm minor, and 

D(A) 



(I + A)~ 1 = I + 



det (J + A) 



Since det (/ + A) is differentiable on the complex Banach space Ji it is analytic 

J-1 „„J -i-l £ „„ :„ 7~1/ A \ 
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Theorem 3.1.1. If A(x) G Ji is holomorphic in a complex variable x, and det(/ + 
A(x)) 7^ for at least one x ^ xq, then (I + A(x))~ 1 is meromorphic in x, with 
poles at the zeroes of det (I + A(x)). 

In the course of the solution of the factorization problem we need to consider 
operators on H± of the following type: 

C gf = Y fk9j-k, j > 0; C~f = fk9j-k, j < 0. 

fc>0 fc<0 

3.1.2. Theorem. For g G G, are Fredholm operators on H±. 

Proof. Let K g denote the operation defined on H by K g f = f 'g, where g G G and 
/ G H. The equation K g f = h can be written as the infinite system of equations 

^2f k gj-k = hj. (3.1.3) 



This system is invertible for g G G since g is bounded and det g(X) = 1 on the unit 
cle, and K g has a bou 
We rewrite (3.1.3) as 



circle, and K g has a bounded inverse on H. 



fk9j-k + J2?k 9 J+ k =h+, j > 0, 

fc>0 fc>0 

J2 fk9-j-k + J2 fk~9k-j =hj j>0 

k>0 k>0 



(3.1.4) 



where f£ = f k , k > 0; f k = /_ fe , k > 0; etc. 

Denoting {/+, /"} as well by /, we introduce the operators 

B 9 f = Yl fk9j-k, Yl fk~9k-j I , T g f = I ^ fk~9j+k, Y fk9-(k+j) I • 

\fe=0 fc>0 / \fc>0 fc=0 / 

By the decomposition (3.1.4), K g = T g + B g , and B g = K g — T g = K g (I — K~ 1 T g ). 
The operators are precisely the components of B g . Since K g is invertible, it 
suffices to show that T g is compact. 

In fact, each of its components is of trace class. Consider for example, the 
operator on H + defined by Tf = T<kfk9k+j- Now an infinite matrix A = \\ajk\\ is 
of trace class if 

Y, i a j fc i < 

j,k>0 

In fact, (cf. [Si]) there is a partial isometry U such that \A\ = U*A, where \A\ is 
the positive definite operator (AA*) 1 / 2 . Since = H^ll? au entries of U* are 

bounded by one in absolute value. Consequently 

||A||i = tr\A\ = ^2u* k a k j < \a jk \. 
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Thus it suffices to show that T,j jk \gj +k \ < +00 where |<7fc| denotes the sum of the 
absolute values of g k . We have 

00 

Yl \ g j+k\ = ^2( n + 1 )\9n\; 

j,k>0 n=0 

and the sum on the right converges since g is analytic in |A| > 1. 

The second component of T g is also of trace class. A similiar argument leads to 
consideration of the series 

J2\9n\(l + \n\). 

n<0 

This sum converges since g is analytic for p < |A| < 1. □ 
Since T g is of trace class, the following is immediate: 

3.1.5. Corollary. The operators are invertible iff det (I — K~ l T g ) 7^ 0. If g 
is holomorphic in x then (C^) _1 are meromorphic in x. 

Let 

G l = {g : g = g(X) e M n (C), det g = 1, g analytic in t < |A| < +00}, 

and define H l , H± and G l _ accordingly. We are now ready to prove: 

3.1.6. Theorem. Let g E G p . The factorization problem 

<P-g = <!>+, g^G, $± e G±, t>p, 

has a unique solution provided det (I — K~^Tg) 7^ 0. Moreover, if g is analytic in 
a parameter x, the factors <p± are meromorphic in x, provided that for at least one 
value ofx, det (I - K^Tg) ^ 0. 

Proof. Writing 

(f>- = I + ^2<t> k X k , </>+ = ^ fc A fc , 

k<0 k>0 

and substituting these expressions into our factorization problem, we get 




For k < we get the system of equations 

9 k + h9k-j =0, k < 0. (3.1.7) 

j<0 

By Corollary 3.1.5 this system has a unique solution in H l _, denoted by 0_, provided 
that det(7 — K~ l T g ) does not vanish. Setting 

<Pk = 9k + ^2 <Pj9k-j, k>0, (3.1.8) 
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we obtain the required factorization provided we can show that <p± belong to G± 
respectively. Since 0_ and g are analytic in |A| > £, so is <j> + = 4>-9- The factor- 
ization shows that there are no negative powers of A in the Laurent expansion for 
0+. Therefore 0+ is entire. 

Note that det0+ = det0_ on |A| = 1 so the sectionally holomorphic function 
given by 

J det</>_, |A| > 1 

\ det0 + , |A| < 1 

is analytic on the extended complex plane. Therefore it is a constant. Since 
0_(oo) = /, it is identically 1, and det <p± = 1. Hence 0+ is entire and 0± G 
Since 0_ does not depend on t it belongs to G l _ for all t > p, hence 0_ G G_ . 

If (7 depends analytically on x, then the K~ 1 T g is an analytic operator valued 
function of x, and the solution 7 of (3.1.7) is meromorphic in x. If g is entire, then 
the poles of <fi- do not cluster in the finite complex x plane. □ 

3.2 The dressing method. Once the factorization of g G G has been obtained, 
a flat connection is reconstructed from the factors by a procedure often referred 
to as the "dressing method" [ZS]. We first summarize the method for the AKNS 
hierarchies with potentials in B since these are somewhat less complicated. 
Given g G G, let m± be the solutions of the factorization problem 

e xzJ+tz k „ g ^ e -xzJ-tz k ^ = m -i m+ _ (3.2.1) 

where J and fi are diagonal matrices. It is immediate from (3.2.1) that 

D(mZ 1 m + ) = z[J 1 mZ 1 m + ] 1 D t (mZ 1 m + ) = z k \\i, mZ 1 m + ] 1 (3.2.2) 

where for the present D = d/dx and D t — d/dt. 
Denote by m the sectionally holomorphic function 



r m+ 


\z\ 


\ m_ 


\z\ 



A simple calculation based on (3.2.2) shows that the sectionally holomorphic func- 
tions 

(Dm)m~ 1 + m(zJ)m~ 1 , (D t m)m~ 1 + m{z k \i)m~ x 

have no jumps across the unit circle, and are consequently entire functions of z. 
Since m_ tends to / as z tends to infinity, their highest order terms in z are 
respectively z J and z k fi; hence, with m = m_, 

(Dm)m~ 1 + m(zJ)m~ 1 =u(x, z) = zJ + q(x, t) (3.2.3a) 
(Dtmjm' 1 + m(z k n)m~ 1 =v k (x, z) = z k [i + B k (3.2.3b) 

where B k denotes a polynomial of degree k — 1 in z. 

Letting ( ) + denote the projection onto the non- negative powers of z, we obtain 

Z J + q = (m-zJmZ 1 )+ = zJ + [mi, J] (3.2.4) 

where mi is the coefficient of z~ x in the expansion of m_, and 



v k (x,z) = (m-Z k ixm_ 1 ) + . 



(3.2.5) 
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3.2.6 Theorem. v k is a universal polynomial in q and its derivatives with zero 
constant term. Specifically, 

k 

3=0 

where the matrix valued functions Fj satisfy the recursion relations 

[J, F J+1 ] = DFj - [q, Fj], F = fi. (3.2.7) 

The connection D — u, D t — v k is flat, and q satisfies the nonlinear partial 
differential equation 

q t = [J,F k+1 ). (3.2.8) 

Solutions of (3.2.8) which are obtained from the factorization (3.2.1) are meromor- 
phic in x and t. 

Likewise, the coefficients of mJ^wT 1 , where me xJx satisifes (2.1.1), are differ- 
ential polynomials in q with zero constant term. 

Proof. We sketch a proof of the first statement based on arguments in [BS2]. The 
wave function m_ is uniquely determined up to right multiplication by a diagonal 
matrix, so F = rri-Z k fimZ 1 is independent of the choice of m_. For fixed xq, 
construct a formal series m_ by requiring the diagonal entries in nij, j > to 
vanish at xq. It is clear that the coefficients of m_ are then uniquely determined 
polynomials in q and its derivatives at x$, so the same is true of the Fj. Since xq is 
arbitrary, the result holds for all x. When q is zero, m_ is necessarily diagonal and 
F = z k fi; hence these polynomials have zero constant term. The same argument 
holds for the Gel'fand-Dikii case, where m_ satisifies m x = [J\,m] + qm and 
F = mJ\m~ x . 

To prove the second statement, rewrite (3.2.3) in the form (D — u)m = m[D — 
zJ), (D t — v)m = m(D t — z k jj). Then 

[D -u,D t - v]m = m[D - zJ, D t - z h fj] = 0. 

Since m is invertible, [D — u, D t — v ] = and the connection is flat. The evolu- 
tion equation (3.2.8) follows in the standard way from this fact and the recursion 
relations (3.2.7). 

Since x and t appear analytically in (3.2.1), the factors m± are meromorphic in 
x and t by Theorem 3.1.6. The meromorphic behavior of q follows from (3.2.4). □ 

The wave function m is said to "dress" the bare connection D — zJ, Dt — z k \i 
to the "dressed connection" D — u, D t — v. 

Flat connections for the Gel'fand-Dikii and associated hierarchies are constructed 
from solutions of the factorization problem (2.3.4), but require a number of alge- 
braic considerations. Namely, one must decide how one is going to split the con- 
stant factors in the Lie group G. By altering the choice of G± one obtains not 
only the Gel'fand-Dikii flows, but their associated flows as well, such as the mod- 
ified Gel'fand-Dikii flows, etc. Moreover, in the case of the Gel'fand-Dikii flows 
themselves, the factor G- in the Riemann-Hilbert splitting is not a group, so the 
dressing method must be modified. This will be done in §5.1. 

The analogs of (3.2.3a,b) are 

(Dm)™' 1 + mJ\m~ x =J\ + q (3.2.9a) 
(Am)m- 1 + mJimT 1 = j\ + B Ujk . (3.2.9b) 
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3.2.10 Corollary. Let q be defined for x G I C K and suppose L has a Baker- 
Akhiezer function for x <E I. Then q has a meromorphic extension to the complex 
plane. 

Proof. For x £ I the wave functions satisfy the factorization problem (2.3.4) (with 
t = 0). Since the left side is entire in x, the factors m± have meromorphic extensions 
to the complex plane, by Theorem 3.1.6. The meromorphic extension of q is then 
given by (3.2.9a). □ 

This result was obtained in [SW] using the r function; the method here applies 
to the AKNS and Drinfeld-Sokolov hierarchies as well. 



4. Algebraic structure 

4.1 Loop groups and algebras. We begin by reviewing the basic ideas of loop 
groups and algebras [K], [PS]. We continue to denote by G the loop group introduced 
in §3.1. Let | • | denote a norm on M n (C). We may assume that the automorphism a 
on diagonal matrices £>, which we introduced in §2.2, leaves | • | invariant: |cr(L>)| = 
\D\. 

From the identity 

xm D , - . i 



it follows that an equivalent representation of G is given by 

oo 

G = {g : 9 = dkJ ^ det 9 = 1 > d k^ V ) 



k= — oo 



where it is assumed that the series converges absolutely and uniformly on S 1 . Given 
g = Tikdk J^, and h = T,khk J\ we have 



9 h = J2 ^ d *AVfc)j J l 



Moreover, G is a Banach Lie group with norm 

Ibll = ^2\dkl 
k 

and Banach Lie algebra 

g = {h : h = ^/ifcJ^; \\h\ \ < +oo, tr =0, for k = mod n}. 

k 

An integer valued grading of g is a direct sum decompostion of g into subspaces 
2k such that [g j7 Qt] C Qj+k- In our case, g has two natural gradings, induced by 
the two bases: 

degi (g k X k ) = k, deg 2 (d k J^) = k. 
The homogeneous subspaces are different for the first and second gradings; let us 
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4.1.1 Lemma. The complex dimension of q\ is n 2 — 1 while 

2 ( n if k 7^ mod n 
dim g k = < 

{ n — 1 otherwise. 

4.2 Lie algebra decompositions. Let 

-l -l 
V ™ = Q 2 j={a:a= a * J t a k E V}. 

j=— n+1 k=—n+l 

It is easily seen that dim V n = n(n — 1) and that 

where U± are the subspaces of strictly upper and lower triangular matrices. 
G and g are contained in the linear space 

B = {a:a = J2 a ^ J t a k E V} = {a : a = ^# fc A fc , g k E M n (C)}. 
fc fc 

This space also has gradings relative to the two bases, as above. We define projec- 
tions Pi and P 2 (denoted respectively by P + and P + in [DS]) onto the subspaces 
of non-negative degree relative to these gradings by 



P 1 a = ^g k \ k , P 2 a = J2 a kJ. 



k 



fc>0 fc>0 

where 

a = ^9k\ k = J2 akJ *- 

k k 

We define a third projection on B as follows. Let D\ be as in Theorem 2.2.10. 

4.2.1 Theorem. There is a projection P3 whose image coincides with that of P\ 
and whose kernel is given by 

ker P 3 = {g : g = ^2e k J^ k , e k E S k } 

k>0 

where 

span [I, ((T-/)- 1 ^!, ... , (a-^-tt-VDx], k < n - 1; 
V k>n. 



Proof. We must show that every a E B can be written uniquely as a = h + v , where 
v E Im Pi and h belongs to the subspace ker P3 defined above. We may assume a 
is homogeneous with respect to the second grading. If deg2 a > 0, then a E Im Pi, 
and a = v. On the other hand, if deg2 a < —n, then a E ker P 3 and a = h. 

That leaves elements in V n . If deg a = —j then a = ajJ^ J for some aj E V. We 
must show that 
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where 

J} = {/:/ = diag(/i,...,/ B ), / 1 = / 2 = ... = / i = 0}. 

It is clear that dim JFj = n — j. We first prove that dim Sj = j and so it suffices 
to prove that Sj fl Tj = {0}. Suppose that 

j-i 

e = u; 7 + ^u; r ((7 - I)~ r D l = 0. 

r=l 

Then 

(a - 7) (j " 1) e = wi((7 - 7) J ' _2 Di + • • • + = 0. 

It is easily checked that (a — I) k D\ E T n -k-i for k < n; and this implies that 
cUj_i = 0. Applying (cr — I)- 7-2 to e we find next that cu_/-2 = and so forth. Hence 
all the coefficients are zero, and dim Sj = j. 

Now let e be as given above and suppose that e E Sj H Tj . We repeat the above 
argument. Applying (cx — 7) J ' -1 to e we obtain 

(a - 7)^- 1 >e = u>i(<7 - 7)J'- 2 Di + • • • + E (a — iy~ x T 3 C J^i. 

Hence u)j-i = 0. Then we apply (a — 7) J ' -2 to e, and find that = and so 

forth. Hence e = 0, and Sj fl Tj = {0}. □ 

The decomposition (4.2.2) can be obtained by multiplying on the right by etk 
with k < j and taking the trace. In that case we obtain the system of equations 

j-i 

tr[e kk a] = u + ^ w r tr[e fcfc ((j - 7)~ r 7>i], l<k<j, (4.2.3) 

r=l 

for the coefficients u>o, . . . ,Wj-i- For future reference we denote the matrix of 
coefficients of (4.2.3) by W. 

The projection P3 is obtained by splitting V n C ker Pi into the direct sum of 
two subspaces. The dimension of ker P 3 fl V n is n(n — l)/2. 

Let us set 

Qj = ker Pj n , 0+ = im P 3 n fl, j = 1, 2, 3. 

Then 

4.2.4 Theorem. 

i) = flj © 0+ /or j = 1,2, 3. 

ii) 0]~, 02~, fli", 82"? an d 03~ are eac ^ ^ e subalgebras of g for all n. 

iii) 03 is a subalgebra of g for n = 2, 3 but not for n > 3. 

Proof. Statement i) follows from the fact that each Pj is a projection. That the 
kernel and range of Pi are subalgebras is immediate from the definition. To prove 
the corresponding fact for P 2 , note first that J\bJ^ 1 = a(b) for any diagonal matrix 
b. Then, for a and b diagonal, 

r„ tJ u rfcl _ /„_7/i,\ i,_fc/„^ d+k 
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It follows that the kernel and range of P 2 are subalgebras, since they are the sub- 
spaces of elements spanned respectively by negative and positive powers of J\. 
To prove iii) it suffices to consider the subspace 



3 HK = sp{e k J x fc , k = 0, . . . , n - 1, e k E £ k } 

For n = 2, this subspace is spanned by J^ 1 . Since [J^ 1 , J^ 1 ] = 0, it is trivial that 
g% is a subalgebra. For n = 3 one must also check the commutator 

[J A - 1 ,(a-/)- 1 D 1 J A - 2 ] = -a- 1 ( J D 1 )J A - 3 - 

The element on the right has degree -3 in the second grading. When n = 3 it 
automatically belongs to q% , since in that case no further constraint is placed on 
the diagonal matrix multiplying J^ 3 - 

For n > 4, however, the condition that it belong to 03 is 

a~ 1 D 1 = al + (3(a - I)~ 1 D 1 + 7 (ct - I)~ 2 D 1 

for some complex constants a, (3, 7. Operating on this equation by a(a — I) 2 , we 
find 

(a 2 (l-/3) + (/3-2- 7 )(7 + l) J D 1 = 0. 

It is not hard to see, from the definition of D±, that these equations are inconsistent. 
Thus g% is not a subalgebra when n > 4. □ 

4.3 Lie group decompositions. The above Lie algebra decompositions have 
their counterparts on the group level. Let 

G+ = {g:geG, g E ImPj}, (4.3.1) 
GJ = {g : g E G, g - I E KerP,}. (4.3.2) 



Then: 



and 



G 1 = {g : detg = 1, ^ = ^^A J , g = 1} 

j<0 

G+ = {g:detg=l, 3 = ^9^}, 

j>0 



G 2 = {g : detg = 1, g = ^2gjJ{, 9j e X>, fifo = / } 

J<0 

= : det^ = 1, g = ^g j J{, gjEV}, 

j>0 



G 3 = {g : detg = 1, g = ^g j J{, gj E £j, g = 1} 

j<0 

G+=G+ 
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4.3.3 Theorem. 

i) G+ = G+. 

ii) G2 is a Lie subgroup of and G x is a Lie subgroup of G 2 . 

iii) G~ and Gj~, j = 1, 2 are Banach Lie groups with Lie algebras q~ and Qj~ . 

iv) G3 is a group only for n = 2, 3. 

This theorem is proved in the same manner as Theorem 4.2.5; the demonstration 
is left to the reader. 

We are now ready to state the main theorem of this section: 

4.3.4 Theorem. For each of the three projections Pj there is a an open dense 
subset Qj C G such that each g G Qj has a unique factorization 

<7 = (0-)-V;, fieGf. 

Proof. For j = 1 this is the well-known Riemann-Hilbert factorization whose proof 
is given in Theorem 3.1.6. Now suppose 3 G fii admits such a factorization with 
^GGp We write 

</>+(A)=4 + (A), j> + (0) = I, v = <j>+(0). 

Let SL°(n) denote the subset of v in SL(n) which have the factorization v = i~ x u 
where t = I + strictly lower and u is an upper triangular matrix. As is well known, 
SL°(n) is a dense open subset of SL(n). When v G SL°(n), g factors as 

where £(f)~ G G^ and ucj) + G 

The domain 2 consists of all g G fii for which v = 4> + (0) G SL°(n). Since fi 2 
is the inverse image of SL°(n) under the continuous mapping g — > v, it is open. To 
show that Q 2 is dense in O l5 consider g G Oi\0 2 - From (3.1.8), 

v = g + J2^9-j- 

j<0 

Let g denote a tangent vector at g to a curve in G. For those curves for which 
g is independent of A, i) = go- The set of such tangent vectors is transversal to 
SL(n)\SL°(n); for example, it contains all strictly triangular matrices. So every 
neighborhood of g contains points in 2 . 

The construction of the third factorization is more involved. We begin by prov- 
ing: 

4.3.5 Lemma. Given a in Gi or G 2 there is a unique constant lower triangular 
matrix t with 1 's on the diagonal such that ia G G% . 

Proof. It suffices to prove the lemma for a G G 2 , since G± is a subgroup of G% ■ 
Let 

00 

a = V"a fc J^ fc , a k G V, a = I. 
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Any constant lower triangular I with l's on the diagonal can be written as 

n-l 

<=Etf- /o = /, 
i=o 

It follows that 

oo I k \ 

fo = E E/i ff_i ( a *-i) j a* 

fc=0 \j=0 / 

where we set /j = for j > n. Thus la e G3 if and only if 

fc 

^^(st-jjefi, fe = 0,...,n-l. 
i=o 

Writing these equations out we get 

/1 + ai = ei 
/2 + + a 2 = e 2 

n-2 

/n-l + E ^" fcCr fc ( a "-fc-l) + a n-l = e n _l. 
fc=l 

But these decompositions are precisely those obtained in (4.2.2). Given a\ the 
first equation may be solved uniquely for fi and e\. Then / 2 and e 2 are uniquely 
determined from a\ and a 2 from the second equation, and so forth. □ 

We now return to the proof of Theorem 4.3.4. We set = Given g e Oi, let 
g = (fii4>i. By the preceding lemma there exists a unique lower triangular matrix 
I with l's on the diagonal such that £^1 £ ^3"- Note that left multiplication by a 
constant matrix of determinant 1 leaves 67 + = 67g~ invariant; and therefore 

</=(WW)» £<pteGf 

is the desired factorization. 

For j = 1,2, the uniqueness of the factorizations is a simple consequence of the 
fact that Gj and 67 + are groups and 67~ fl 67 + = {/}. For j = 3, however, 63 is 
not a group for n > 3, and we have to argue differently. 

Let us assume therefore that g = ((^>3~) _1 ^> + = (03 ) _1 (/> + . Let £, respectively I 
denote the values of (/)% and at infinity; I and i are lower triangular matrices 
with l's on the diagonal. Writing (p^ = £v$ and ^3 = £v 3 , we have 

g = v^ii- 1 ^) = v^Cr 1 ^), (4.3.6) 

where now 

v 3 , v 3 eg^, rVt r^GG+ 

The factorization (4.3.6) of g is of type 1, and is therefore unique, by what we have 
already proved, so v 3 = v 3 . By Lemma 4.3.5 there is a unique matrix I such that 

such that hio = £11* <= G~ Thus P = £ \ (h~ = d>7 . and rh+ = thf . I~l 
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5. Nonlinear Flows 

5.1 The Gel'fand-Dikii flows. In the preceeding section we introduced three 
different factorizations of elements of G. These three factorizations lead to the 
Gel'fand-Dikii, modified and "potential" Gel'fand-Dikii flows. In this section we 
discuss the construction of these different flows from the different factorizations. 
We begin with the Gel'fand-Dikii flows themselves. 

5.1.1 Theorem. Let g belong to G and let m± be the solutions of the factorization 
problem 

e ^ J ^ tJ 1)g(\) e -^ J ^ tJ 1) =mZ 1 (x,t,X)m + (x,t,X) (5.1.2) 
where m± G Gf. Then 

(Dm)!!!, -1 + mJ\m~ x = J\ + q, m = m±, (5.1.3) 

where q = q(x,t) G B satisfies the Gel'fand-Dikii equations for (/c,n) 

Proof. By the results of §3 and Theorem 4.3.4 the factors m± of (5.1.2) are mero- 
morphic functions of x for fixed t and of t for fixed x. By the arguments of §3.2, 
(Dm)m -1 + mJ^m -1 , where m is the sectionally holomorphic function m±, is an 
entire function of A. 
Since m_ G G% , 

m- = S ^ j m j J~\ mjESj. 

j>0 

In particular, m_ ~ / + £, as A — > oo, where I is strictly lower triangular. Since t 
depends on x, (Dm-)mZ 1 is bounded as A — > oo while 

m-JxmZ 1 — J\ =[m_, JxlmZ 1 

J2( m J ~ a ( m i)) J l~ J m ~- 



vi>2 



m_ 



The above expression is bounded as A — > oo, and therefore so is 

(Dm-)mZ 1 + m-JxmZ 1 — J\- 

It follows by Liouville's theorem that this expression is independent of A. 

Let us call it q(x,t). From the expansion for m_ we conclude that q may be 
expanded as 

n-l 

q = ^QkJZ k , 

fc=0 
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Writing (5.1.3) in the form Dm + mJ\ = (J\ + q)m and expanding q and m 
in powers of J\, we obtain a sequence of recursion relations for q k (noting that 
mo = I): 

k 

Dm k + (I - a)(m k+ i) = ^qjO~ J (m k -j) = qom k H h<?fc, (5.1.4) 

3=0 

for < k < n - 1. 

Since m_ G G 3 , mi = cuio-T, and <?o = — c)(toi) = 0. From (5.1.4) we get for 
fc = 1 

qi = mi + (I - cr)(m 2 ) 

=Dlu 10 I - (a - I)(u2qI + u>2i(a - I)~ 1 D 1 
=Du 10 I - u) 2 iDi. 

Since qi G T\ its first entry vanishes, so W21 = -Du^io, and 

qi = nDu) 10 e nn . 

We now prove, by induction, that in fact qj is a scalar multiple of e nn for all 
1 < < 77. — 1. Assume this is true for all integers less than j. Since the m s are 
diagonal matrices, qimj-i + • • • + ?j_imi on the right side of (5.1.4) is a scalar 
multiple of e nn , so we need to prove that 

Dnij + (I - a)(m j+ i) = dje nn (5.1.5) 

for some constant dj. Since qj G J-j the right side of (5.1.4) belongs to Tj, for 
k = j, and so 

tr e 88 [Drrij + (I - a)(m j+1 )] = 0, s<j. (5.1.6) 

Writing 

j-i 

mj = u j0 I + ^2u jk (a - I)~ k D u 

k=i 
j 

TOj+i = Wj + i >0 / + ^CU i+ i )fc ((T - I)~ k D u 
k=l 

we find, after some calculations, that (5.1.6) reduces to the linear system of equa- 
tions 

j-i 

tr e ss (Du j0 I - Uj+x^Di) + y~](Du} jk - u j+ljk+1 )tr e ss (a - I)~ k Dx = 0, s < j. 

k=i 

This is the homogeneous system of equations corresponding to (4.2.3); and as we 
saw in the proof of Theorem 4.2.1, that system has a unique solution. Therefore 
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and (5.1.5) is proved. 

We now prove that q satisfies the (n, k) Gel'fand-Dikii equations. Since G% is not 
a group, we cannot use the Lie algebra decomposition that worked for the AKNS 
hierarchy. Instead, we argue as follows. The solution of the factorization problem 
satisfies the simultaneous equations (3.2.9a,b). Therefore we have 

{D-J x - g)V = 0, (D t - G fc )V = 0, (5.1.7) 

where Gk = J^+-B n ,A: is the right side of (3.2.9b), and tp = mex.p{xJ\+tJ^}. From 
the first equation we see that ifi is a Wronskian: ip = \\D 3 Vi-\ \\, i, j = 1, . . . , n, 
and that the entries of the first row of ip satisfy Lv = Xv. If xq is a regular point of q, 
then these equations are also satisfied by ifj = ip(x, X)tfj~ 1 (xo, A), and ^(xq, A) = I. 
We shall henceforth assume this is the case and drop the tilde. 
By the second equation each Vj satisfies 

where the G\j denote the entries of the first row of Gj~- Since Lvj = Xvj, we 
can replace X s v by L s v. The result is that each vj satisfies an equation of the 
form vj = PkVj where Pk = D k + . . . is a differential operator of order k and is 
independent of A. The simultaneous equations Lvj = Xvj, and vj = PkVj imply 
that (L — [Pk, L])vj = 0; and moreover, D m Vj(xo, z) = bj m for < j, m < n — 1. 
Let 

k 

[Pk,L]=Y j r j {*)D j - 
3=0 

Then 

k 

= (t - [P k , L]j Vj(x , z) = Uj(x ) - rj(x ) - r m (x )D m v j (x , z), 

m=n 

where the sum over m > n is vacuous if k < n. Since derivatives of v of order 
n or greater can be replaced by A, the summation is a polynomial in A of degree 
at least one. Since the above equation holds identically in A, we conclude that 
the summation vanishes identically, [Pk, L] is of order less than n, and L — [Pk, L] 
vanishes at xq. Since xq was arbitrary, it vanishes at all regular points of q. 

This implies that [Pk, L] is of order n — 2. It follows (cf. Proposition 2.3 [DS]) 
that the coefficients of Pk are differential polynomials of the coefficients of L and 
that Pk is a linear combination of the operators L J ^ n for j < k: 

k 

Pk= £ Cj (L^) + , 

for some constants Cj. 

To show that Pk is precisely L+ n , it suffices to prove that all the Cj, j < k, 
are zero. We do this by showing that for j < k the coefficient of D J in P k is a 
differential polynomial in the entries of q with zero constant term. The coefficients 
of Pk come from the first row of Gk, and this matrix is in turn given by (3.2.9b). 
From the lower triangular structure of the leading term of m G G^ , the first row of 
Gk is identical to that of Pi(mJ^m _1 ). But the entries of mJ\m~ x are differential 
polynomials in q with zero constant terms by Theorem 3.2.6. □ 
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5.2 The Modified GePfand-Dikii flows. The so-called modified Gel'fand-Dikii 
flows are a special case of the AKNS hierarchy, and are defined as follows. Writing 
L = (D + u n )(D + u n -i) •••(£) + ui), the equation Lv = Xv can be written as a 
first order system 

(D - J x + q)i> = 

where q is the diagonal matrix with entries u±, . . . ,u n . ( tr q = since the coefficient 
of D n ~ x vanishes.) 

Recall (§2.1) that A~ 1 J A A Z = zJ a , where J a is diagonal. Since q is diagonal, 
A~ x qA z = q is independent of z. Hence we can write the isospectral problem as 

(D - zJ a + q)$ = 0. 

This is the isospectral problem for an AKNS system, whose flows we constructed 
in §3.2. 

Working in the original basis, we obtain the (n, k) modified Gel'fand-Dikii flow 
as [D x , D t ] = 0, where 

F = mJ\m~ x , and [•]+ denotes the projection P%, i.e. the projection onto non- 
negative powers of J\. 

The modified Gel'fand Dikii flows are obtained from the second factorization: 

5.2.1 Theorem. Let g e G and let m± be the holomorphic solutions to the fac- 
torization problem 

et { x j x+ tj k x ) ge -i{xj x +u k x ) = m= i m+5 m± e G ± 

Then 

q(x, t) = (Dra)m -1 + mJ\m~ x — J\ 

is a diagonal matrix of trace zero and satisfies the modified Gel'fand-Dikii equations 
for (k,n). 

Proof. We repeat the first part of the proof of Theorem 5.1.1. Again we may 
conclude that 

(Dm)m~ 1 + mJxm' 1 - J x (5.2.2) 

is an entire function of A, where m is the sectionally meromorphic function defined 
by m+, m_ in the regions |A| < p, |A| > p respectively. This time 

ui- = S ^iTLjJ\'' , mo = I, tTij G V. (5.2.3) 

Again m_ = I + £ + 0(X~ 1 ) as A — * oo and so the expression in (5.2.2) is bounded 
as A tends to infinity. As before, we denote it by q, where q depends only on x and 
t. Inserting (5.2.3) into (5.2.2), we obtain q as the constant term in 

771- J\Wr} 

= (I + m 1 J- 1 + .. .)J X (I - mi J" 1 + ...)- J x 
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where the dots denote negative powers of J\. Thus 

q = (/-(t)N GD . 

To show that q satisfies the modified Gel'fand-Dikii (n, k) flow, we consider the 
connection D — J\ — q, Dt — vt, where Vk is given by the right side of (3.2.9b). As 
in the proof of Theorem 3.2.6, this connection is flat. This time G~^ is a group, so 
we obtain the Lie-algebra decomposition 

v k =P 2 [(Am)m -1 + mJ^m" 1 ] 
=P 2 (mj£m- 1 ) 



k-j 



J=0 

By conjugating by the constant matrix A a everything is expanded in powers of zJ a , 
(Fj is diagonal, so it commutes with d(z)) and we see that this is precisely the k th 
AKNS flow. □ 

5.3 The "potential" Gel'fand-Dikii flows. Factorizations in Gf lead to a third 
structure for the potential. This time m G G± and so has the expansion 



m 



E m i t 
AT' mo = L 



3=0 



The arguments of the previous section are repeated, this time using the projection 
Pi. From (5.2.2) we obtain 

Pi((Dm)m~ 1 + mJxm' 1 ) = P^mJxm -1 ) = J\ + [mi,e„i] 

and one finds that 



/ TOi„ 

m 2n 








rrinn-mn -m 12 



° \ 



-m ln J 



where rrijk are the entries of the coefficient matrix m\. The dressing argument 
goes through as in the previous section. The following is related to the general 
theory of Drinfeld and Sokolov. As they show, all reductions of the n th order 
equation Lv = Xv to a first order system are related by a gauge transformation, 
and the entries of L are differential polynomials in the coefficients the potential. 



5.3.1 Theorem. The entries of q are differential polynomials of precisely n — 1 
independent functions. 
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5.3.2 Lemma. For m G G l we have the representation 

mo = I, nij G V (5.3.3) 
where the diagonal matrices mi, . . . ,m n satisfy the additional constraints 



m = m j J 

3=0 



j 

X ' 



trmjekk = j = 1, . . . , n - 1; k > j + 1 
trrrij — for j = mod n. 



(5.3.4) 
(5.3.5) 



Proof. From the form of J A 1 , namely 



Jx 1 



/0 ... 

1 

1 

V 



o / 



we see that the constraints (5.3.4) are necessary and sufficient that the constant 
term in m be the identity matrix. The constraint tr m n = is a consequence of 
the condition det m = 1 . □ 



Now write q in the form 



n-1 

( l = ^2 ( lj J x 3 

J=0 



where 



q = u (en - e nn ), qj = Uj(ejj - e nn ) + Vj(ejj + e nn ), q n -i = v n -ie nn . 

Substituting this form of q, together with the representation (5.3.3) into the equa- 
tion Dm = [J\, m] + qm, we obtain the recursion relations 

mi - a{mi) = q , 
Dmi + m 2 - cr(m 2 ) = qomi + qi, 



n-1 



Drrij + m j+ i - a(m j+ i) = ^ Qk<* k (m j - k ). 

k=0 

We solve these equations recursively. From the first we get mi = wo e n an d 
Qo = Mo(en — e nn ). Taking the trace of the second equation we find tr qi = 2v\ = 
Duq — Uq. With this choice of vi the diagonal matrix m 2 is uniquely determined 
by the second equation and the constraint (5.3.4). This process continues up to 
j = n—2. At each stage the parameters Vj are determined in terms of the preceeding 
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The equation at j = n — 1 is 



Dm n -i +m n - a(m n ) = qom n -i H h qVz-i- 



(5.3.5) 



The first n — 1 rows of (5.3.5) plus the trace condition tr m n = uniquely determine 
the entries of m n as differential polynomials in Uj, Vj for j = 0, . . . , n — 2 The last 
row of (5.3.5) then determines i> n -i 

We thus see that the entries of q are differential polynomials of the functions 

u , ...,M n -i- □ 
For n = 2,3 we get 



This class of potentials leads to a third class of flows, the "potential Gel'fand- 
Dikii" flows. These flows may be obtained as a reduction of the potential Kadomtsev- 
Petviashvili hierarchy [Sz]. 
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